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Abstract 



The gravitational strength of the central singularity in spherically symmet- 
ric space-times is investigated. Necessary conditions for the singularity to 
be gravitationally weak are derived and it is shown that these are violated 
in a wide variety of circumstances. These conditions allow conclusions to be 
drawn about the nature of the singularity without having to integrate the 
geodesic equations. In particular, any geodesic with a non-zero amount of 
angular momentum which impinges on the singularity terminates in a strong 
curvature singularity. 
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I. INTRODUCTION 



The theorems of Hawking, Penrose and others predict the occurrence of space-time sin- 
gularities in a variety of interesting physical situations The singularities which have 
received the most attention over the last years are those which occur in gravitational col- 
lapse and the initial cosmological singularity. It seems fair to say that our understanding 
of these singularities remains at a preliminary stage; little is known about generic 4-d col- 
lapse and correspondingly, generic inhomogeneous cosmological singularities. However much 
progress has been made on the understanding of these singularities under certain simplifying 
assumptions, e.g. the assumption of spherical symmetry for black holes or of homogeneity 
for cosmological singularities. See for example the reviews of and 0]. 

It is in this context that we analyse a particular feature of singularities, namely their 
gravitational strength [^], under the simplifying assumption of spherical symmetry. This 
continues the work initiated in [^, but here we concentrate on the central singularity. The 
notion of the gravitational strength of a singularity was first introduced by Ellis and Schmidt 
0] with the aim of distinguishing between singularities which destroy objects impinging upon 
them and those which do not. A formal mathematical definition was given by Tipler, based 
on the familiar idea of modelling an object's physical extension using Jacobi fields along its 
world-line [Q. Thus let 7 : [to,0) — >• M (with tangent /c°) be an incomplete causal geodesic 
running into a singularity as the parameter t — 0^. We define the sets of Jacobi fields 
Jt^_Tio < ti < as follows: Jt^ = : (i), (ii), (iii)}, where 

(tt) r(ti) = 0; 

{ill) D^C + Rbcd'k^k'^C = 0. 

(Covariant differentiation along 7 is represented by D.) Note that the elements of Jt^ are 
space-like vector fields, and so we may refer to their norm, ||^|| = {gab^°'^''Y^'^- Given any 
three (for time-like geodesies; two for null) independent elements of J^^, a volume element 
along 7 is constructed by taking the exterior product of the corresponding 1-forms, and a 
volume V by taking the norm of this 3-form. Then the singularity is said to be gravitationally 
strong if for all such volumes V, we have 

liminf y(t) = 0. 

According to Tipler's definition, the singularity is said to be gravitationally weak if this 
condition does not hold. Thus any object with world-line 7 will inevitably be crushed by a 
gravitationally strong singularity. To emphasize the role played by 7 here, we will refer to 
geodesies terminating in strong or weak singularities. 

As pointed out recently and independently by Nolan and by Ori [0, this definition 
ignores some singularities which would destroy objects impinging upon them and so needs 
a brief addendum. Firstly, the definition of a strong singularity ignores the case where 
V diverges to infinity in the approach to the singularity. Subject to the strong energy 
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condition and Einstein's equation, is a convex function of t, and so cannot diverge in a 
finite amount of parameter time [§. However, there are situations where the strong energy 
condition is violated while the weak and dominant energy conditions are satisfied. In such 
a case, convexity of V is not guaranteed and so one should allow for the possibility of V 
diverging. Secondly, as pointed out by Tipler, the volume form may stretch infinitely in one 
direction while shrinking to zero in another in such a way that its norm V remains finite 
overall. Such a 'spaghettifying' effect clearly signals the end of an observer's history. (An 
observer falling radially into the singularity at the centre of a Schwarzschild black hole suffers 
infinite stretching in the radial direction and infinite crushing in the tangential directions. 
The net effect on his volume V is that it is crushed to zero. An explicit example of a 
situation where the radial stretching and tangential crushing are exactly cancelled when one 
calculates V was given in 0.) Such situations should also be included in the definition of 
strong singularities. This may be done in a logical and succinct fashion following Ori [0: a 
singularity is said to be deformationally strong if it is either (i) Tipler strong (i.e. strong in 
the sense of the paragraphs preceding this one) or (ii) if for every ti, there exists an element 
of Jt^ which has infinite norm in the limit t 0" . A singularity is said to be deformationally 
weak if it is not deformationally strong. 

Our aim is to give useful geometric conditions for the occurrence of deformationally strong 
singularities. We focus on a particular class of singularities; those which occur at the centre 
of spherically symmetric space-times. Numerous analyses have predicted the occurrence 
of such singularities inside spherical black holes PHT5[|, as a consequence of gravitational 
collapse |[l6|-p0| and in cosmological models We exploit the available symmetry to 



develop a condition which is necessary for a singularity to be deformationally weak. In 
conjunction with an existing necessary condition for a singularity to be Tipler- weak [Q, 
we find very severe restrictions on the existence of weak central singularities. In fact it 
will be shown that any non-radial causal geodesic which approaches r = terminates in a 
deformationally strong curvature singularity. For radial causal geodesies and in particular 
situations (e.g. assuming a particular matter distribution or a space-like singularity), we can 
show that these restrictions are violated, i.e. the singularity is deformationally strong. In 
most cases, we present our necessary conditions as inequalities which must be satisfied along 
geodesies running into the singularity. However, in many cases we can draw conclusions 
without having to integrate the geodesic equations. Thus strong curvature singularities can 
be predicted at the level of the cuvature tensor rather than the geodesies themselves. 

In the following section, we study causal geodesies and the volume V along them in 
spherically symmetric space-times. This allows us to give the result referred to above for non- 
radial geodesies and to present our main result on radial geodesies in the form of necessary 
conditions for a central singularity to be deformationally weak. Applications are then given 
and further comments are given in a concluding section. We emphasize throughout the use 
of invariant quantities. 



II. GEODESICS IN SPHERICAL SYMMETRY 

We write the line element in double null form; 

ds^ = -2e"2^ dudv + r^dQ^, (2.1) 
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where / = f{u, v), dfl"^ is the hne element of the unit 2-sphere and r = r{u, v) is the radius 
function of the space-time (which is a geometric invariant). u,v are null coordinates and 
the form is invariant up to v — >■ Vi{v), m — >• Mi(-u). A singularity will be referred to as 
central if it occurs at r = 0. In the coordinates of ( p.l|) , the Ricci tensor has non- vanishing 
components 

Ruu = — 2r (r u„ + Ivuf.u) (2.2a) 
i?,, = -2r-i(r,, + 2r,/,„) (2.2b) 
Ruv ~ 2r {t ^uv f ,uv) (2.2c) 
Ree = csc^ OR^^ = 1 + 2e2^(r,„r,, + rr,„,). (2.2d) 

Insofar as it is possible, we will describe any curvature tensor terms that we encounter using 
the following invariants; e^^ Ruu, R^, the Misner-Sharp energy p2 | 



E = ^(l + 2e2^r„r,), 

the Newman- Penrose Weyl tensor Coulomb component (calculated on a principal null tetrad) 

e^/ E 
*2 = — (r ™ + rf^uv) - 

and the Ricci scalar R. A useful feature of E is that it offers a simple description of the 
trapped and untrapped regions of a spherically symmetric space-time; the point x G M 
lies on a trapped (untrapped, marginally trapped) 2-sphere iff x = 1 — 2Er~^ is negative 
(positive, zero) at x ||23|] . 

Given an arbitrary geodesic 7 in spherical symmetry, the coordinates of the 2-sphere 
{6, (p) may be chosen such that the motion proceeds in the hypersurface 9 = 7r/2. Thus the 
tangent to an arbitrary causal geodesic may be written as 

k = u— + v— + Lr-^ — , 2.3 

OU OV 0(p 

where we have included the conservation of angular momentum, = L = constant. The 

overdot indicates differentiation with respect to the parameter t. The remaining geodesic 
equations are 

- 2e-^^ iiv + = -e, (2.4a) 

u - 2/>2 + ^2g2/ ^~3^^^ ^ Ah) 

V - 2f^y + L'^e^f r-^r^u = 0, (2.4c) 

where e = -|-1 for time-like geodesies and e = for null geodesies. 

In a previous paper p, we studied radial causal geodesies and were able to obtain a 
useful decomposition of the volume V. We found that V = \axy\, where a is the norm of a 
radial Jacobi field and x, y are the norms of two mutually orthogonal tangential Jacobi fields 
along the geodesic. The key to obtaining this decomposition is the fact that any three such 
Jacobi fields (which vanish at t = ti) provide a basis for J^^. This complete decomposition 
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of V is not available in the general non-radial {L ^ 0) case, but a useful part of it is. This 
partial decomposition relies on the following facts. 

As we will see, there is always a Jacobi field along 7 of the form 

^=^^"'§0- (2-5) 

x{t) must satisfy a certain ODE along 7; see below. In the time-like case, the other two 
elements of a basis for Jt^ must have non-zero components in the 2-space orthogonal to 
both k and ^, and indeed may be taken to lie in this 2-space. We consider two such basis 
elements, fj and (. Now introduce an orthonormal tetrad, parallel propagated along 7. An 
arbitrary tangent vector / orthogonal to k will only have components on the three spatial 
vectors of the tetrad; let these components be a = 1,2,3. For three such vectors, the 
corresponding 1-forms satisfy 

||Ii AI2 AI3II = det 

where the columns of the matrix are the tetrad components of the given vectors. For the 
case where the /(j) are elements of Jt^, this matrix is a constant matrix multiple (which we 
will call a transition matrix T) of the matrix 

r = 

By orthogonahty, the volume associated with risl^ = ||^||(||^7AC||). In the case of relevance 
to us, where the /(,) are independent elements of J^^ , the transition matrix is non-singular 
and so we obtain for the V of relevance, 

V{t) = det{T)Vr. 

The key point here is that |x| = ||^|| appears as a factor of V(t), and so if x(t) is degenerate 
(i.e. X — or 00) in the limit as the singularity is approached (t — 0~), then the singularity 
must be deformationally strong. A similar argument holds in the null case. 

We turn to the derivation of the equation satisfied by x{t) in ( |2.5| ). It is easily verified 
that the vector field r~^d/dd is parallel propagated along 7 and has unit norm. Thus ( p.5|) 
satisfies the geodesic deviation equation iff 

x5^e = -Rw/k^'k'^x. (2.6) 

We find that 

RbefkV = (-(2r-V,„/,„ + r-V^^Jn^ - (2r-V,,/,, + r-'r,,,)v^ 

+ sin^ ^(1 + 2e^^ r,„r - 2r-^r^uvUv) 

This term is controlled solely by r{t). The evolution of r along 7 is given by r = r^uii + 
r^yi). The second derivative can be worked out and simplified with the use of the geodesic 
equations, resulting in 
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Comparing the last two equations along 7 (on which 9 = vr/2), we see that the equation 
( p.6|) becomes 

^+(^-^)^ = 0- (2-7) 

The existence of the angular momentum term plays a vital role as we will see in the following 
section. We treat the radial (L = 0) and non-radial (L 7^ 0) cases separately, dealing with 
the latter first. 



^A^"^) + z^y = o- (3-1) 



III. NON-RADIAL GEODESICS 

We may assume without loss of generality that L > 0. Defining y = r~^x, the equation 
may be written in the self adjoint form 

dV dt' ' r2 

We can obtain the asymptotic behaviour of x in the limit as the singularity is approached 
as follows. First, we move the singularity out to infinite parameter value. The origin 
and temporal orientation of the affine parameter /proper time t has been fixed so that the 
singularity at r = is approached as t — > 0~. We define s = —t~^, so that the singularity is 
approached as s +00. Defining R{s) = r(t), where r(t) means r{u, v)\u=u(t),v=v{t), i-e. this 
indicates the dependence of r on the parameter t in the solution of the geodesic equations, 
equation (|3.1| ) becomes 

(s'RVy + = 0' (3-2) 

where the prime denotes differentiation with respect to s. Thus the equation is of the form 

{p{s)y')' - q{s)y = 0. (3.3) 

We can use the Liouville- Green asymptotic formula to obtain the leading order behaviour of 
y and hence x as s ^ 00. We quote in full the following theorem which appears as Theorem 



2.2.1 of H 



Theorem 1 Let p and q be nowhere zero and have locally absolutely continuous first deriva- 
tives in an interval fa, 00). Let 



pq \ P , 

and let 



00) (3.4) 



{p-^/\-^/\pqy)' eL{a, 00) (3.5) 

Let Re((^-|-w2)^/2) have one sign in [a, 00) where w = {pq)' / {^pq)- Then (\3.'^ has solutions 
1/1,1/2 with asymptotic behaviour 

1/1,2 ~ (M)-'/'exp £(^ + w^)^/^d^ , {s ^ 00). (3.6) 
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We now turn to the application of this theorem to the equation We have here 

2 d2 ^ 

P=sR, ^=-^r 

Recall that we are working under the hypothesis that R{s), satisfying the geodesic equation, 
approaching i? = as s ^ oo. Thus there exists a real a such that neither p nor q are zero 
in the interval [a, oo). We have 

p' = 2sR^ + 2s^R! = 2r^r^ + 2r, 

q' = 2L\s-^R-^ + s-^R^^R!) = 2L\t\-^ + t\-^r). 

The requirement that these functions be locally absolutely continuous on [a, oo), i.e. that 
they be absolutely continuous on compact subsets of [a, oo), is very weak. This would follow 
from local boundedness of p', q' which itself will follow from the assumption of the existence 
of the geodesic on the interval, i.e. from the fact that f is defined and bounded away from 
the singularity. Also, 

pq = -L^, 

which is constant, and so conditions ( p.4|) and (|3.5| ) are automatically satisfied. Note also 
that w = 0. Thus the hypotheses of the theorem are satisfied and the conclusion ( |3.(j| ) yields 
(on removing a constant unimodular factor) 

Writing this in terms of the affine parameter /proper time t, we obtain for the two real 
independent solutions Xi^2 of (^TTI) 

xi(t)~L-VV(t)cos (/^rff), 
x,(t)~L-VMt)sin^'-A_rft'^. 

Now let J"'{t) = x(t)r^^6g be a Jacobi field which vanishes at some arbitrary time ti < 0, 
i.e. x(ti) = 0. Since both of the independent solutions Xi^2 of ( p.7|) approach zero as t — > 0^, 
we conclude that the particular linear combination which gives the present x{t) will also 
approach zero in this limit. Thus we see that in every case, x{t) — »• in the limit as the 
singularity is approached. This proves the following result, which we note is independent of 
any energy conditions. 

Proposition 1 Let j be a non-radial causal geodesic in a spherically symmetric space-time 
{M,g). If 'J runs into the centre r = in finite parameter time, either in the past or in the 
future, then 7 terminates in a deformationally strong curvature singularity. 

7 



(3.7) 
(3.8) 



Thus any singularity, naked or covered, which is reached by a non-radial causal geodesic 
is deformationally strong. These geodesies have not been widely studied and deserve some 
attention. It would be of interest to know, for example, if those space-times which admit 
radial geodesies with past endpoints on a central singularity - i.e. which admit naked 
singularities - also admit non-radial geodesies with the same, or indeed if there are non- 
radial geodesies whose futures terminate at the past endpoints of the naked singularity 
geodesies. If such geodesies exist, then these singularities should be considered genuine; 
their existence has a destructive effect on certain observers in the space-time. We hope to 
address the question of the existence of such geodesies in future work. A useful starting 
point would be the invariant equation for the evolution of r along the geodesic, i.e. equation 
(|]3D below. 



IV. RADIAL GEODESICS 

In the radial case, the angular momentum term L vanishes and ( p.7|) reads 



rx — rx = 0. (4.1) 

The unique solution (modulo an irrelevant constant factor) of this equation satisfying x{ti) = 
is 

x{t)=r{t) (4.2) 

Notice then that if the singularity is non-central, a;(0~) is non-zero and finite. For a central 
singularity, we have the following 



Lemma 1 Let 



I 



ds 



ti r2(s) 



(i) If the integral I converges, then x ^ as t ^ . 
(a) If I diverges, then 

lim xit) = lim — (t). 

This is a straight application of I'Hopital's rule. From this and the comments above, we 
have the following useful corollaries. 

Corollary 1 If a radial causal geodesic 7 terminates in a deformationally weak central 
singularity, then along 7, lim^^o- '^(^) is non-zero and finite. 

Corollary 2 Let the conditions of corollary |i| be satisfied. Then there exists cq > such 
that 

r(t) ~ Co|t| as t — s> 0^. 
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This yields a useful necessary condition for deformational weakness of the singularities 
under consideration. We use corollary |l] above and some established results on weakness of 
singularities |p to derive a new necessary condition for deformational weakness of central 
singularities reached by radial geodesies. This condition seems unlikely to be satisfied in 
many circumstances. We require the following preliminary basic results. 

Lemma 2 Let a E C{0, b] for some b > 0. Suppose that 

a{s)ds 

converges for all r G (0, e)(6 > e > 0). Then limr-,o ra{r) = 0. 
Proof: Let r be fixed and define (3 on [0,r] by 

rr 

j3{x) = f a{s) ds. 







Then j3 G C[0,r] n C^(0,r]. By Taylor's theorem, for every x G [0,r], there is an ri G (x, r) 
such that 

= /?(r) +/?'(ri)(x-r) 
= (r — x)a{r-i). 

In particular, 

Hir) := r a(s)ds = p(0) = m(r,) 
Jo 

for some r^, G (0, r). By hypothesis, H{r) exists and is finite for all sufficiently small r. Thus 

= lim|i7(r)| 

1 — >o 

= lim |ra;(r^,)| 
> hm 

= lim |r*a(r*)| > 0. 

Replacing r^, by r in the last line gives the required result. 
Lemma 3 Let r(t) satisfy the differential equation 

r = a{r) 

on (0,6] where a G C{0,b] for some b > 0. Suppose that \imt-,or{t) = and that lim^^o^ 
exists and is finite. Then limr^o ra{r) = 0. 

Proof: For sufficiently small r and for t > 0, we may integrate the differential equation to 
obtain 
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f2(t) = 2 r a{r')dr' + r'^{to), 

Jro 

where r(to) = ''"o ^ ^- Then by hypotheses, 

/ a{s) ds = hm / a(r') dr' 

J ro r— >0 J J.Q 

exists and is finite. Now apply lemma ^ to obtain the result. 

Next, we recall a result of Clarke and Krolak (a direct consequence of their corollary 
3 For this we note that if a singularity is deformationally weak, then it is Tipler- 

weak. Note that this result assumes the strong energy condition and Einstein's equation (or 
equivalently, the time-like and null convergence conditions). 

Lemma 4 If a causal geodesic 7 terminates in a deformationally weak singularity, then 
along 7, 

lim t^i?44 = 0, 

where -R44 := Rabk°'k^- 



In order to use lemma 0, we note the following. We have r = r uii + r yV. Differentiating 



again, using the geodesic equations (p.4|) and grouping terms appropriately, we get the fol- 



lowing form for f (for generality, we include angular momentum in the following expression): 

r''r = -ii?44 + e(4 + 2^2 - ^) + L'r'\\ + S^^)- (4.3) 
2 r-^ 12 

The right hand side is to be viewed as a function of the parameter t (i.e. it is assumed that 
the geodesic equations are solved), which is smooth for < \t\ < 6 for some 6 and singular 
at t = 0. (The degree of smoothness is not particularly significant; continuity is sufficient. 
However there is very little restriction in assuming a higher degree of differentiability for 
\t\ > 0; we are interested in singularities occurring at t = 0.) In the present situation, 
corollary ^ applies and so we can use the inverse function theorem to write the right hand 
side of ( [4. 3D as a function of r which is continuous on (0, b] for some positive b. Thus lemma 
H applies. 

Combining lemma ^ with corollaries ^ and 0, we obtain our main result. 

Proposition 2 If a radial causal geodesic 7 terminates in a deformationally weak central 
singularity, then along 7 

limr^i?44 = lim rf = 0. 
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The usefulness of the result comes from the fact that on the one hand we have two 
independent conditions which must be satisfied by weak central singularities, and on the 
other, the parameter t does not appear explicitly in the relevant quantities. These features 
are emphasized by using (|4.3| ) and the following equation for R^^: 

■ E i?, 9 9/E i? , 

- + *2 - - + - + ^2 - TTT, 

r"* r-^ 24 



Ru = RuuU^ + RvvV^ + 2e(- + *2 - -) + 4L2r-2(- + ^2 - — ). (4.4) 



An immediate consequence of equations ( |4.4| ), ( [4. 3D and proposition^ is the following: 



Corollary Z If a radial causal geodesic 7 terminates in a deformationally weak central 
singularity, then along 7 

limer2(4 + 2^2-^)=0. (4.5) 
r^o 12 

Notice that this result is vacuous for radial null geodesies; this is a consequence of the fact 
that these are principal null directions in spherical symmetry. On the other hand, we see that 
the strength of the singularity approached by a radial null geodesic is completely controlled 
by the behaviour of i?44 in the limit as the singularity is approached. Indeed we can give 
the following result. 

Corollary 4 A radial null geodesic which reaches r = in finite parameter time terminates 
in a deformationally weak singularity if and only if there exists e > such that rR^^ is 
integrable as a function of r on [0,e]. 

Proof: In the radial case, the volume element along the geodesic is a constant multiple of 
x^, where x satisfies ( [4.2|) . This follows from the fact that the second Jacobi field may be 
chosen to be y(t)r^^ esc 96'^, with y also satisfying ( |4.2| ). See 0. Thus by Lemma 1, the 
singularity is deformationally weak if and only if r is finite in the limit as the singularity is 
approached. But using (^73|), we can write 



- rr'R44{r')dr' + f\to), 

J rn 



where tq = r(to)- The result follows immediately from this integral. 

This result shows that the converse of Lemma 4 (with t replaced by r) is true for radial 
null geodesies. This may be useful as this lemma has been used widely in studies of radial 
null geodesies emanating from central singularities; see e.g. |1^JT^. We now have the useful 



converse, that if the condition used to demonstrate deformational strength of the singularity 
fails, then the singularity is necessarily deformationally weak. 

In the following section, we analyze the conditions r'^R^^ and rf — subject to var- 
ious assumptions. In several cases, we show how deformationally strong central singularities 
may be identified without having to integrate the geodesic equations. 
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V. APPLICATIONS 



In this section, we assume the following situation obtains: there exists a radial time-like 
geodesic 7 which runs into r = in a finite amount of parameter time. Then the origin of 
the parameter t along the geodesic may be translated so that r(0) = 0. We set aside the 
issue of the existence of such geodesies. 

We assume that the dominant and strong energy conditions are satisfied by the energy- 
momentum tensor of the space-time and that Einstein's equation holds. The dominant 
energy condition states that T"feZ^ is past-directed and causal for any future-directed time- 
like 1°". Then in particular Tabl°'l^ > for all causal 1°". We choose units so that Einstein's 
equation is Gab = ^T^Tab- Given this equation, the strong energy condition is equivalent 
to Rabl°'f > for all causal In particular, the dominant energy condition yields the 
following inequalities: 



E ^ R 
^ ■ ■ 12 



■.= ^ + ^2 + ^>0. 



The approach we take here is to derive general results about the strength of the singularity 
without integrating either field equations or geodesic equations. We will try to derive general 
results based on certain geometric assumptions, and otherwise, restrict to particular matter 
models. In this case, we will focus on two important and widely studied scalar field 

in a source free electric field, and a perfect fiuid. In the former case, the Ricci tensor is given 
by 

i?,6 = 2V,0V,0+^E„,, (5.1) 

where is the scalar field, Q is the constant electric charge, Eab = —Qah + 2r^Sab and Sab is 
the standard metric on the unit 2-sphere. We have R = 2g°'^'Va4'^b4' ^^id 

We note that the dominant and strong energy conditions are automatically satisfied by this 
matter distribution. THe neutral case can be studied by setting Q = 0. Burko has 
shown that under the assumption of spatial homogeneity, the central singularity in this 
model is deformationally strong. 

For a perfect fiuid with fiow vector u"', energy density p and pressure p, we have 

{Sny^Rab = (p + p)uaUb + ^(P - P)9ab, (5.3) 

i?= 87r(p-3j9), (5.4) 

E A-nn , , 

- + vl/2 = -^. 5.5 

r-^ 6 

The dominant energy condition requires p + p>0, p — p>0. 

We now proceed to investigate the consequences of proposition |^ and corollary ^ under 
the following cases. Unless otherwise stated, we make no assumptions about the mat- 
ter distribution, other than that the dominant and strong energy conditions are satisfied. 
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Throughout the remainder of this section, asymptotic relations, hmiting values etc. refer to 
the limit as r — along a geodesic which terminates in a deformationally weak singularity. 

For radial time-like geodesies, L = 0, e = 1 and so the conclusion of corollary |^ is that, 
along a radial geodesic terminating in a deformationally weak central singularity. 



-'(f + 2^2-^)-0. (5.6) 



For convenience, we give 



i?44 = RuuU^ + RvvV"^ + 2(^ + ^2 



E ^ R. 
r-^ 6 



1 1 „ , E R . 

r-'r = --Ru + i-,+2^,--). 

It is difficult to make any general statements without making further assumptions. How- 
ever ( |5.(j| ) may prove to be a useful condition to use to check the strength of certain singu- 
larities. We can make some progress in the case where 

4 + *2-f>0. (5.7) 

Then in addition to (|5.6| ), we must have 

r, , E ^ R. 

We can write these as (any two of) 

r\z^ - -) ^ 0, r\z^ - 4) ^ 0, r\z^ + 8^2) ^ 0. 

If any one of < or i? < or \I'2 > holds, then these give 

r^^a, —,r'^R 
r 

as conditions which must hold at a weak singularity. Notice then that the singularity must 
be untrapped (1 — 2E /r — 1 > as the singularity is approached). That is: 

Corollary 5 If a radial time-like geodesic 7 terminates in a deformationally weak central 
singularity and if the inequality (\5. 1\) and at least one of the inequalities 



E <0, R<0, *2 > 
hold in the limit as the singularity is approached, then the singularity is untrapped. 



It is worthwhile investigating the consequences of (|5.7| ) for the two matter models men- 
tioned above. For a scalar field in an electric field, this is equivalent to 
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which may be described as the case where the electric field dominates the scalar field. Then 
we conclude that 



7 = ^ + ^(1)' 

Notice that the singularity must be naked. This should be considered in conjunction with 



the result of [1^, where evidence is presented that in this case, the singularity must be time- 
like but that the evolution tends to avoid this situation (i.e. the scalar field dominates). For 
the case of a neutral scalar field ((5 = 0), ( ^.7|) becomes 



Setting Q = in the previous trio of equations, we see that in this case the singularity is 
untrapped. 

For a perfect fiuid ( p.7|) reads 

P > 0, 

and for a deformationally weak central singularity subject to this condition, we conclude 
that 



r^p^O, 27rp = — + o(r-2). (5.8) 



Ori and Piran considered self-similar collapse of a perfect fiuid with a barotropic 
equation of state (which must necessarily be of the form p = 7p, 7 constant). Every such 
space-time includes a central singularity at r = 0,t = where t is an orthogonal time 
coordinate, fixed by demanding that it measures proper time of an observer at the regular 
centre. This point is referred to as the origin. It is found that r^p = D{x) where x = r/t 
is the similarity variable. Thus by ( |5.8|) (which applies if < k < 1/3), a radial time-like 
geodesic running into the origin terminates in a deformationally strong singularity provided 
limD(a;) 7^ along the geodesic. The existence of such geodesies is readily demonstrated 



using the results of [0. The corresponding result for outgoing radial null geodesies was 
proven by Ori and Piran; this includes the interesting case of future-pointing outgoing 
radial null geodesies originating at the singularity i.e. the case of a naked singularity. The 
present result shows that these singularities will destroy an observer impinging upon them. 
Deformationally strong spacelike singularities have also been detected in the gravitational 
collapse of a perfect fiuid (with and without an electric field) under the assumption of spatial 
quasi- homogeneity p5| . 

Recall from above that the present case (|5.7|) includes dust. The analysis which follows 
gives a non-trivial example of the demonstration of the deformational strength of a singu- 
larity which does not rely on solving the geodesic equations. Thus we investigate the latter 
condition in (WM, i.e. 
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r\27vp - ^) ^ 0. 



For the case of marginally bound dust, we have 

E' 



Attp 



where the prime denotes differentiation with respect to a spatial coordinate x which labels 
points in the slices orthogonal to the fluid flow. The Einstein equations in this case yield 
E = E{x) and 

r%x,T) = ^EiToix)-Ty, 

where the fluid flow vector is Ua = —'VaT. See [0,0 for details. There is freedom in the 
choice of the coordinate x which allows x — » X(x). This may be utilised by taking r = x on 
the initial slice r = 0. This choice has the advantage of specifying to{x): 

ro(x) = 3^x3/^. 

The central singularity in this model appears at r = to(x); that singularity for which we 
also have x = is of interest for studies of cosmic censorship. Jhingan and Joshi argue that 
the appropriate form for E is 



E{x) = Y.En 

n=0 



where the En are constants. {E is not required to be analytic and this is not supposed to 
be implied by the form above. In what follows, all we require is that E ~ Eqx^^x and 
that this relation is differentiable at x = 0. This form for E is the most general that ensures 
a finite, non-singular initial state for the matter.) With this assumption, a straightforward 
calculation gives 

^ r^' E'r{To~T) + 2T[fEr' 

This quantity definitely diverges at a central singularity r = (tq — r) = 0,x 7^ 0; such a 
singularity must be deformationally strong. In the case r = (tq — r) = x = 0, we have 

SEox' + ^E'J'E^' 



r (Airp - 2—) - ,„ 

r-^ Q5/30-1/3 Z7'l/'^/'^_ -T-^5/3 _ O-I/6Q-I/3 Z?^'/o : 



35/32-1/3^1/3^^^ _ ^)5/3 _ 2-1/63-1/3^17/6^-1^^^^ _ ^^2/3 

Again, this quantity (generically) diverges, giving a strong curvature singularity. We em- 
phasise that it was not necessary to integrate the geodesic equations in order to reach this 
conclusion. 



15 



VI. CONCLUSIONS 



We have exploited the symmetry properties of causal geodesies in spherically symmetric 
space-times to study the effect of a central singularity on an observer who impinges upon 
it. We have been able to demonstrate the destructive effect (deformational strength) of the 
singularity in several cases, and have given very finely tuned conditions which must hold in 
order that this destruction need not occur. In particular, any non-radial geodesic approach- 
ing a central singularity must terminate in a deformationally strong central singularity. This 
may be understood as follows. The focussing term i?44 will include the term R^^cl?, which 
by the conservation of angular momentum equals L?R^^r~'^. This introduces strong curva- 
ture along the geodesic which contributes to the destructive effect. Of course this does not 
include the vacuum case, and so we conclude that the angular momentum must also cause 
a significant amount of shear to develop along the geodesic which, via the Raychaudhuri 
equation, contributes to the strong focussing effect. 

The main advantage of our approach was that it did not require the integration of the 
geodesic equations; it was possible to predict the deformational strength (or weakness) of 
certain singularities by calculating the Riemann tensor rather than its tetrad components 
P). Of course one needs to address the issue of the existence of geodesies which run into 
the singularity (i.e. the question of the existence of the singularity) in the situations studied 
above. However this can often be done without having to obtain the detailed and subtle 
information required to apply the results of |g] (see, for example, |^ for a thorough appli- 



cation of these results to the null weak Cauchy horizon singularity in spherical black holes). 
For example, if = is satisfied along a causal geodesic - as is the case for a neutral scalar 
field - then ( |4.3| ) reads 

rr = -"^{RuuU^ + Rvvv"^) - e— + ^(1 - 3—). 

In a trapped region, —E/r < —1/2, and so f < 0. Assuming the absence of singularities 
away from r = 0, this is sufficient to ensure that the geodesic runs into the centre. If L 7^ 0, 
this will be a deformationally strong singularity. 

The analysis here was made possible by the assumption of spherical symmetry. One 
would expect similar results in space-times with hyperbolic and plane symmetry. It may also 
be possible to extend the applicability of the idea of determining the nature of singularities 
from simple geometric quantities to more general situations, e.g. axially symmetric space- 
times or homogeneous cosmologies. We note that in this vein, significant progress has 
been made recently on the issue of the connection between a well-behaved metric and weak 
singularities . 
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